January 20, 2013 8:46 WSPC/INSTRUCTION FILE ijmpa-HOTDl 



International Journal of Modern Physics A 
© World Scientific Publishing Company 



o 

(N 
o 

o 



Perturbative Hamiltonian constraints for higher order theories. 



S. A. Martinez, R. Montemayor 
Instituto Balseiro and CAB, Universidad Nacional de Cuyo and CNEA, 8400 Bariloche, 

Argentina 



Oh. 



> 
m 

o 



X 



L. F. Urrutia 

Instituto de Ciencias Nucleares, Universidad Nacional Autonoma de Mexico, A. Postal 70-t 

04510 Mexico D.F., Mexico 

Received Day Month Year 
Revised Day Month Year 

We present an alternative method for constructing a consistent perturbative low energy 
canonical formalism for higher order time-derivative theories, which consists in appliying 
the standard Dirac method to the first order version of the higher order Lagrangian, 
augmented by additional perturbative Hamiltonian constraints. The method is purely 
algebraic, provides the dynamical formulation directly in phase space and can be used 
in singular theories without the need of initially fixing the gauge. We apply it to two 
paradigmatic examples: the Pais-Uhlenbeck oscillator and the Bernard-Duncan scalar 
field with self-interaction. We also compare the results, both at the classical and quantum 
level, with the ones corresponding to a direct perturbative construction applied to the 
exact higher order theory, after incorporating the projection to the space of physical 
modes. This comparison highligths the soundness of the present formalism. 
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1. Introduction 

Higher order time-derivative (HOTD) theories, including nonlocal theories, have a 
long history in physics. They continually reappear in new models of physical inter- 
est, particularly in the form of effective Lagrangians that describe small corrections 
to well established theories. For example, nonlocal effective field theories emerge 
when high energy degrees of freedom are integrated oulP. HOTD terms also appear 
in higher de rivativ e gravitj^, nonconmutative field theorjl^, models derived form 
string theorji^EEl^ effective models for meson nucleon interactioni^, etc. Our con- 
cern in this work is the discussion of HOTD as perturbations of well established 
standard second order theories (to be called precursors), instead of considering them 
as fundamental ones. Then it is appropriate to view the resulting HOTD description 
as an effective field theory, valid within a given energy range, where a perturbative 
description is perfectly acceptable, besides of constituting a natural setting for ex- 
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tracting physical predictions. 

Classical HOTD theories introduce more degrees of freedom than their precur- 
sors, which is better reflected in their Hamiltonian structure. One characteristic of 
these theories is the Ostrogradsky instability, the existence of unphysical runaway 
solutions, not expandible in powers of the parameter codifying the HOTD contribu- 
tion. At the quantum level the changes in the canonical structure produce important 
differences, associated to several pathologies. Even when the HOTD terms are con- 
sidered as small corrections to their precursors, as it is done here, their effect is 
qualitatively significant. We understand here the number of degrees of freedom as 
the number of initial conditions that must be given to fully determine the behavior 
of the system at any time. Naturally, this depends on the degree of the time deriva- 
tive in the equation of motion, irrespective of the fact that we might, for example, 
have only one coordinate x{t). 

The most traditional approach for a canonical formalism for HOTD theories 
was developed by Ostrogradsky in 1850 in the context of regular theories. This 
construction is equivalent to finding a second order Lagrangian by introducing the 
appropriate auxiliary variables in the HOTD theories via Lagrange multipliers, and 
by applying the Dirac methoci^. The canonical formalism thus obtained highlights 
the well known problems of these theories, such as instability, Hamiltonians un- 
bounded from below and lack of unitarity at the quantum level. The Ostrogradsky 
method was extended to HOTD singular theories by NesterenkJI^. 

However, the HOTD terms usually appear in the Lagrangian as small cor- 
rections, labeled by a small parameter 7. For this reason, perturbative schemes 
have been developed which identify the correct low-energy degrees of freedom 
contained in the the ory and avoid the high-energy ones that produce all the 
inconvenienc^^HmH A typical recent example of this is the Myers-Pospelov 
modeff^, whereby HOTD operators of dimension five are introduced as small cor- 
rections to standard electrodynamics in order to describe possible minute signals 
of Lorcntz invariance violation. An analogous situation occurs in Lorentz violating 
extensions for the Standard Modeff^J^ which generalizes the previous situation to 
the standard model of particles plus gravity. HOTD theories are also present in 
the analysis of precision tests of electroweak interactions Those research topics 
are the subject of a considerable number of experimental observations with ever 
increasing precision. A very useful toy model for explorii ig t he different aspects 
of such perturbative methods is the Bernard-Duncan fielcP^, a g eneralization to 
field theory of the Pais-Uhlenbeck quantum mechanical modell-^, a paradigmatic 
example of HOTD theories. 

There are two main approaches for reducing a higher order Lagrangian to one 
containing only the low energy mode s. One of t hem is based on the use of field 
transformations containing 

derivativePESEOEI]. 

The other approach is based on 
a Lagrangian which contains HOTD up to a certain order. The dynamics can be 
reduced to a second order one by introducing perturbative constraints. Within a 
Lagrangian perspective, these constraints become naturally generated by equations 
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of motion. One of the explored possibilities is to eliminate the higher order deriva- 
tives in the Lagrangian using the equations of motion. Given that in general the 
equations of motion can not be introduced in the Lagrangian without distorting the 
variational principle, this approach can be applied only to certain special cases to 

00 0^ 

o btain an approximate second order Lagrangian - Another possibility, explored 
in^I^, is to project the Lagrangian constraints into the phase space of momenta and 
coordinates, and consider these projections as Hamiltonian constraints, restrict- 
ing the Ostrogradsky Hamiltonian. The key this approach is that the perturbative 
Lagrangians constraints are truly pro tectable and form a set of second class con- 
straints. But, as has been shown ini^, Lagrangian and Hamiltonian constraints are 
not directly related, and in general not projectable in the case of gauge theories. If 
there is a gauge symmetry this approach is not directly applicable, because the first 
class constraints are not projectable. For this reason it is necessary to first fix the 
gauge at the Lagrangian level to proceed with the construction. 

Once the role of the Lagrangian constraints in the reduction of the phase space 
is appreciated, the emphasis of the problem shifts to finding more efficient methods 
of calculating the iterative steps which are required to obtain the sought approx- 
imation to a given order in 7. One of these alternatives is proposed in Refs. |4j 
m where the construction is performed basically in the coordinates-velocities (CV) 
space, avoiding an explicit projection of constraints. The perturbative Lagrangian 
constraints are directly implemented on the Noether energy of the HOTD theory, 
which become the time evolution generator in the constrained CV space. To find 
the canonical structure this constrained energy is considered as the corresponding 
Hamiltonian, and the dynamics is written in terms of generalized brackets among 
coordinates and velocities. Imposing that these brackets reproduce the perturbative 
equations of motion to the order considered, the corresponding algebra is deter- 
mined. This step requires the solution of a set of second order differential equa- 
tions, which become very involved beyond the first order in 7 and whose solution 
involves a good amount of guess work. The algebra in the CV space is subsequently 
rewritten in terms of canonical coordinates and momentum. Ref. [26| includes a very 
clear review of this work, together with that concerning the general problem of 
HOTD systems. There, the need to add constraints in order to make perturbative 
sense of HOTD theories is also emphasized. Constrained HOTD theories become 
free from the diseases that plague unconstrained HOTD ones. This has also been 
previously remarked in Ref. [271 Finally we mention Ref. [TU where the use of iter- 
ative solutions of the equation of motion is focused on the obtention of the final 
symplectic form in terms of the variables q, q, or equivalently ttq in the notation 
of Ref. im The substitution of equations of motion into the Lagrangian, which is 
normally forbidden, is justified here in virtue of the detailed construction of Ref. 
I12| which shows that this substitution really amounts to strongly imposing a set of 
second class constraints which lead to Dirac brackets. In particular the method pro- 
posed in Ref. [T2]can be applied only to regular higher order Lagrangians. The above 
methods are based on the projection of the Lagrangian constraints into the phase 
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space of momenta and coordinates, and subsequently considering these projections 
as Hamiltonian constraints. But, as has been shown in Ref. 24, Lagrangian and 
Hamihonian constraints are not directly related, and, what is more, in general they 
are not projectable. In principle, this restricts the applicability of this approach. At 
least, in the case of higher order gauge theories it is necessary to implement a gauge 
fixing before applying these approaches. 

Previous arguments made clear that HOTD theories, in particular when con- 
sidered as corrections to standard ones, require the imposition of perturbative con- 
straints at a given level . This implementation is closely related to the basic question 
of which are the appropriate Feynman rules to calculate a given low-energy process. 
To this end it is necessary to understand the difference between HOTD theories 
in the context of effective field theories and their use as HOTD theories per se. In 
Ref. [18] such a difference is explored in the context of the Bernard-Duncan modeff^. 
In brief, given the full HOTD Lagrangian, the Mathews' theorem leads to the 
Feynman rules as read directly from this Lagrangian, which will imply the use of 
the full propagator as is proved in Ref. [161 Nevertheless, as we know, this theory has 
all the problems of the HOTD theories. Hence, to make it consistent as a pertur- 
bation of the standard scalar field, we should expand the propagator in powers of 
the coupling constant associated to the higher order term. With this manipulation 
we obtain a perturbative expression for the propagator of the usual scalar field, and 
because of the disappearance of the high-energy poles, the ghost degrees of freedom 
are no more present. The important point to be stressed here is that a perturbative 
expansion of the propagator of the higher order theory provides us with a reference 
to test the soundness of an effective Lagrangian. 

In this paper we propose an alternative con stru ction, which amounts to the 
application of the well established Dirac method^, to the first order version of 
the HOTD theory, augmented by additional perturbative Hamiltonian constraints. 
This construction leads to the perturbative canonical formalism in a very system- 
atic and simple way, working from the beginning in the corresponding phase space. 
The construction is purely algebraic and does not require the solution of any sys- 
tem of differential equations. Also, the method can be directly applied to gauge 
theories without the need of initially fixing the gauge, as would be the case in the 
constructions of Refs. (H [121 

The organization of the paper is the following. In the next section we intro- 
duce the concept of perturbative Hamiltonian constraints, which allows the appli- 
cation of the Dirac method to obtain a consistent canonical formulation, exact to a 
given order in the perturbative parameter. There we also show that the additional 
(perturbative) contraints can be consistently implemented in the Hamiltonian for- 
mulation. The method is illustrated in the third section with the construction of a 
canonical formalism of arbitrary order in the perturbative parameter for the Pais- 
Uhlembeck oscillator. The fourth section applies the method to the Bernard-Duncan 
theory "'^^ with a </j^ interaction, both at the classical and quantum level. Finally, 
the fifth section discusses the two particle quantum scattering in the model of the 
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previous section in order to compare the exact canonical theory, which contains 
ghosts, with the well defined perturbative one obtained with our method. We show 
that our canonical formalism recovers the results obtained by using perturbative 
propagators and ruling out the unphysical states in the exact, pathological, theory. 
The last section contains a summary of the main results and some general com- 
ments. In the Appendix we consider a simple model to illustrate how our method 
compares with the proposal in Ref. 21 

2. Extended Dirac approach: Hamiltonian perturbative constraints 

Any HOTD Lagrangian can be rewritten as a first order one by introducing an 
adequate number of auxiliary variables, via Lagrange multipliers. Once this is done 
we can calculate the corresponding momenta using the usual definition. Some of 
these relations allow us to write a restricted set of velocities as functions of the 
coordinates and the momenta, while others yield constraints. The usual procedure 
to consistently define the dynamics is the use of the Dirac approach, and in this 
way we obtain the canonical formalism, which in the case of a HOTD theory is 
plagued with several pathologies. In our case, the HOTD terms turn out to be 
scaled by a parameter 7, which we assume small, in such a way that the HOTD 
contributions are considered as perturbations over a standard theory. Here it is 
possible to go directly to a perturbative Hamiltonian construction, because some 
of the relations that come from the definition of the momenta and the constraints 
are inhomogeneous in the perturbative parameter. These relations, multiplied by 
a power of 7, provide perturbative constraints valid up to a given power of 7. 
These new constraints can be treated as Hamiltonian ones, and added to the set of 
original primary constraints. Thus, besides those primary constrains generated by 
the definition of the momenta, this approach requires the introduction of additional 
primary constraints further imposed according to the order of 7 to which we decide 
to incorporate the HOTD corrections. From now on the procedure follows as in the 
usual Dirac analysis of constrained systems. 

In the following we show that this is a consistent way of considering such 
constraints. To make this point clear, let us assume that we start with a La- 
grangian L — L{q, q), where p = dL/dq, and that we want to introduce an external 
phase-space constraint, i.e. one not generated by the definition of the momenta, 
f (QtP) = 0- The most straightforward way of doing this is by reformulating the 
theory in an enlarged space {q,p), where the auxiliary variable p corresponds to 
the momentum. If the second order Lagrangian is regular, the procedure for this 
enlargement is given in Ref. I29( and leads to a well known first order Lagrangian of 
the form 

L=p,q'^Hiq,p), (1) 
where H{q,p) is the corresponding Hamiltonian. Here we can impose the additional 
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constraint by using a Lagrange multiplier. Thus we set 

L^p.ct -H[q,p) + \f{q,p). (2) 

Next we apply the Dirac procedure to this Lagrangian, with , pi and A considered 
as coordinates in the extended space. The definition of the momenta yields three 
primary constraints 

T^qi -Pt^O, TTp ~ 0, TTA ~ 0, (3) 

with the extended Hamiltonian 

Hp^ H {q,p) - \f{q,p) + (tt^^ - p^) + Viir^ + wtta, (4) 

where u^^VijW are arbitrary functions. The consistency under time evolution of the 
primary constraints fixes two of the arbitrary functions 

dH ^ df 

oqi oqi 



{TTgi ~ -J> = --^ + A— , (5) 



K,i/}.O^V.|f-Af, (6) 

where ~ indicates a weak equation. These equations generate a secondary constraint 

{^A,i/}~0^/~0. (7) 
Thus, at this level we have the Hamiltonian 

Hn ^ H iq,p)-Xfiq,pH[^^ - A^) - p.)- - Af ) (8) 

with the set of constraints 

TT,, -K~0, TT^^O, ^A~0, /~0. (9) 

The Poisson bracket of an arbitrary function M(q,p) with Hd can be simply 
written as 

{M,HD}^{M,Hy ~X{Mjy , (10) 

where {,}' is the Poisson bracket in the {q,p) subspace. Thus, the consistency 
condition for the secondary constraint / is 

{f,HD}^{f,Hy . (11) 

If {/, iJ}' ~ there are no more constraints, and hence nx is a first class 
constraint. It generates an orbit of equivalent configurations and we can choose 
any point of this orbit by imposing a gauge fixing. It is convenient to take A = 0. 
Besides this, the first two constraints in Eq. ^ are second class. We finally get, in 
the reduced space defined by the gauge fixing and these two second class constraint, 
that the dynamics is described by the usual Hamiltonian Hj^ = H {q,p), with only 
the constraint we want to impose, f{q,p) — 0. 
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If {/, iJ}' is not weakly zero there is a new constraint 
{{f,HD}.HD} = {{f,H}\HD] = {{f.H}\H]' -\{{f,H}' , /}', (12) 
whose consistency condition fixes the remaining arbitrary function w 

[{{f,H}' ,H]' -X{{f,H}' ,f]' ,Hd] = 

{ {{/, H)' , H}' ,h}'-X { {{/, H}' , H}' , /}' - A {{{/, H}' , /}' , h} 
+ X' {{{/, H}' , /}' , /}' -w{{f, H}' J}' . (13) 
Now we have the following set of second class constraints 

TTg,-p,C^O, TT^C^O, TTA-O, / ~ 0, (14) 

{{/,iJ}',i7}'-A{{/,i7}',/}'^0. (15) 

Using the first three constraints in and the last one in to partially reduce 
the phase space we get 

with / ~ 0. In both cases, {/, i?}' = or {/, i?}' ^ 0, the final result corresponds 
to considering the original Hamiltonian in the {q,p) phase space, plus the external 
constraint as a Hamiltonian one in the Dirac approach. 

H{q,p) plus f{q,p)^0 O Hmrac^ H {q,p) +u f {q,p) . (17) 

This construction shows that in the regular case the inclusion of external con- 
straints involving momenta is equivalent to considering these constraints as primary 
ones in the Dirac formalism. The key for this demonstration is the construction of 
a first order Lagrangian with a {q,p) configuration space. In the case of a singular 
Lagrangian, this first order Lagrangian with a {q,p) configuration spac e can also 
be constructed following an extension of the scheme given in Lanczos^, developed 
in Ref. 1301 Using this construction and following the preceding discussion, we can 
show that in any case external constraints involving coordinates and momenta can 
be incorporated as primary Hamiltonian constraints in the framework of the Dirac 
method. In particular, this justifies the insertion of the perturbative constraints in 
the Dirac formalism, together with the ones naturally generated by the definition 
of the momenta. 

In the following we illustrate the proposal by applying it to two systems: (i) the 
Pais-Uhlenbeck oscillator and (ii) the higher order scalar theory discussed in Ref. 
1161 plus a self-interacting term. 
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3. The Pais-Uhlenbeck oscillator 

We implement the proposal in the framework of the Pais-Uhlenbeck oscillator, which 
is a regular system. For arbitrary theories described by a first-order Lagrangian we 
have a well-defined consistent canonical approach, given by the Dirac method. For 
this reason, our first step will be to rewrite any higher order theory in terms of a 
first order Lagrangian by introducing auxiliary degrees of freedom together with 
the corresponding Lagrange multipliers. In this way we obtain a first order singu- 
lar Lagrangian to which we can apply the Dirac method to obtain the canonical 
formalism. When we eliminate the auxiliary variables in this canonical formalism, 
by implementing the corresponding second class Hamiltonian constraints, we re- 
cover the standard Ostrogradsky approach EUHl^ An alternative approach in the 
case of the Pais-Uhlenbeck oscillator based on complex canonical transformations 
plus subsequent reality conditions is presented in Ref. [31] 

The second-order Lagrangian defining the Pais-Uhlenbeck oscillator is 

i-^-^-;*^. (18) 

2 2 2 ^ ' 

The introduction of the additional degree of freedom z, via the corresponding con- 
straint, leads to the first-order Lagrangian 

i=Y-^-^7i^ + A(z~i), (19) 

with coordinates x, z, X. The corresponding equations of motion are 

z -I- 7z + A = 0, (20) 
-uj^x + X = 0, (21) 
z-i = 0. (22) 

They imply only one perturbative primary constraint at order 7" 

7"(z + A)=0, (23) 
which leads to the secondary constraint 

7" (i + a) ^ 7" (i + Lu^x) = 0. (24) 
One can easily verify that eliminating z, and A from the resulting equations of 



motion produces the equation of motion for x directly obtained from Eq. (jl8l) . Now 
that we have a first-order theory, we can construct the canonical formalism using 
the Dirac approach. The corresponding momenta pk : p, tt, k are 

dL 

p= — = -X^ (72 = X+pc^O, (25) 
ox 

TT = — = -7z, (26) 
dL 

K=^ = 0^cri=K~0. (27) 

dX 
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This set of equations gives the transformation of the momenta into coordinates and 
velocities. The equation (|26p defines a perturbative constraint of order n 



$" = 7"7r ~ 0, (28) 

which means that we take 7™ =0, m > n. In this way the velocity i is expressed 
in terms of the associated momentum. There are two primary constraints, and one 
perturbative primary constraint is imposed. The primary Hamihonian is 

? ? 99 

TT Z LO X 

H ^ \ Xz + u{p + X) + vj"tt + WK. (29) 

27 2 2 

The time evolution of the exact primary constraints and the perturbative one yields 

{p + X,H} = -ij?x + w ^ w = u?x, (30) 

{k, -ff} = z — u — )■ u = z, (31) 

7" H} = 7" (z + A) ^ = 7" (z + A) ~ 7" (z - p) ~ 0. (32) 

Two arbitrary functions are fixed, and a secondary constraint is generated. The 
Hamiltonian becomes 

2 2 22 

77 = - ^ + + zip + W7"7r + w^xK. (33) 

27 2 2 

The consistency conditions for the secondary constraints generates a pair of towers 
of constraints of decreasing order in 7, starting from the initial ^I^" 

^~0^ ^^^0^ .... (34) 

at at at 

according to 

{H, $"-™} = 5-"-™, {H, ^r"-"'} = (35) 

The general exp ress ion for the constraints can be written in terms of the Fi- 
bonacci polynomials 

[(™+i)/2] . -u 1 _ A 

Fn.iy)^ J2 (-1)' ™ i jy^ ^o(y) = l, F^,{y)^l, (36) 

where [a] means the integer part of a, and are 

$" = 7'V, (37) 
= 7"^™ (^^^m iW) - pF,n-i (W)) < m < n, (38) 

*"""=7""™(7r7^m-l(7^')-^7^'^m-2(7w')) 1 < TO < n. (39) 

The Fibonacci polynomials satisfy 

Fm+i{y) = Fmiy)-yFr^^i{y), Fo{y) ^ 1, F_i(y) = 1, (40) 
which imply the recurrence relations 

^f"-™ ^ + a;2^"-'"+i < TO < n, ^--1 = ^-"+1 = 0, (41) 

= ^$"-"'-1 + ,j2^n-m+i 1 < m < n, = $"+1 = 0, (42) 
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for the constraints. 

Using these recurrence relations we can express the complete set of constraints 
as proportional to the last ones, ^I^^ and Thus we finally get 

with 

$0 = (7^2) - X 7w2f„_2 (7^') . (46) 

It is important to emphasize that the final perturbative constraints <I>° imply 
all the remaining ones up to 4*", and for this they are taken as the independent 
perturbative constraints to the order considered. 

For simplicity, the illustration of the complete Hamiltonian construction is re- 
stricted here to the case n = 1, defined by the external perturbative constraint 
77r ~ 0. The non-perturbative constraints (|25l27p remain valid for all orders. In this 
case the two independent perturbative constraints become 

$° = Xl = TT - juj'^x, (47) 
= X2 = z- {1 + W)P- (48) 

The system is now second class, and it is very easy to compute from Eq. (j33p the 
Dirac Hamiltonian in the constrained space 

H=^p^ + ^uj^l-joj^)x^, (49) 
together with the corresponding Dirac bracket 

{p, x}^ = !-(!+ 7^') (1 + 7c^') 7^' - (1 - 7c^') • (50) 
It is clear that x and p are not canonical conjugate variables. Instead we can use 

X ^ X, p — > p = (l — 7a;2) p, (51) 

which satisfy 

{p,x}^ = {l-juj'){p,x}^c,l. (52) 

Thus we get 

H=lp^ + loj^l+uj^j)x^ (53) 

to first order in 7. Now we can read the frequency of the oscillator in the Hamilto- 
nian, which coincides with that obtained from the perturbative equation of motion 
to the order considered. 



(43) 
(44) 
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4. Higher order scalar field theory: Bernard-Duncan field with a 
ip'^ interaction 

We now consider a higher order field theory, given by a Lagrangian density which 
is a self-interacting generalization of the Pais-Uhlenbeck oscillator 

£ = -i^n^-i^mV + |'pnV-f/, D^a^s^. (54) 

The HOTD equation of motion is 

□V3 + m^(^ - 7n^(y9 + a^s^ = 0. (55) 

Using an iterative procedure for the Pais-Uhlenbeck model, the above equation 
reduces to the following second order equation 

□<y9 + (l - 7TO^) Lp+{l- A"iw?) aip^ - i^a^-p^ + 6ja(f {df,(fd''(f) = 0, (56) 

to first order in 7. 

We deal with this model following the prescription in the previous subsection, 
so that the Lagrangian ([51)1 is first rewritten as a first-order one with respect to 
the time derivatives, by introducing an auxiliary variable ip = dof with a Lagrange 
multiplier A 

£' = iv^2 _ 1^ (^2 _ A) ^ + 7^AV' + |v,A2^ -^^^ + x(^^-^) + (57) 
The definition of the canonical momenta establishes 

TT^ = 7V', (58) 
and gives two Hamiltonian primary constraints 

XI = TTA ~ 0, (59) 
X2 = - A ~ 0, (60) 

In the first place we will consider this theory without any approximation in 7, and 
in the next subsection it will be reformulated using the Hamiltonian constraint 
approach. 

4.1. The -f-non-perturbative formalism 

Here we construct the canonical formalism using the Dirac approach. The primary 
Hamiltonian density becomes 

H = (1 - 27A) (m^ - A - 7A2) Lp+jLp'^+\ij}+w (vr^ - A)+u7rA, 

^ (61) 
where w and u are arbitrary functions. The consistency conditions for the primary 
constraints only fix two of the arbitrary functions 

{vr^ -\, H}^- (m^ - A - 7A2) - - u ~ 0, (62) 
{7rA,i?}--V + ^«-0, (63) 
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where H = J dP'xH. No further constraints are generated, and thus the Hamiltonian 
density becomes 

^ " ^""^ - V' + 7A^ V + (I - tta) (m^ - a - 7A2) ^+'^^p^ + - aTrx^p^ 

(64) 

The primary constraints are second class 

{Xi{x),X2ix')}=S{x-x'), (65) 
and therefore we can directly use the reduced Hamiltonian density 

n = + - ^V' (1 + 27A) V + (m^ - A - 7A2) ^ + (66) 

with the Dirac brackets 

{if {x, t) , TT^ (y, t)}jj = {ip {x, t) , TT^ (y, t)} = 5 (a; - y) , (67) 
{tp {x, t) , TT^, {y, t)}jj = {tp {x, t) , TT^ (y, i)} = 5 (a; - y) . (68) 

The theory discussed by Bernard and Duncan corresponds to a = 0. We will consider 
now this case, which can be treated in an exact way. To do this, it is useful to 
decompose the field in modes with well defined frequency and with covariant 
normalization 

^p{x,t) = J ^y|= (flpc'^P-''-'"''*) + a],e-^(P-''-'"''*)) 



+ ^== ( fopc'^P-^-^''''^ + foTe-^^P-^-'^''*^ ) I . (69) 
The characteristic frequencies result 

1/2 



2 ^1 + 47771,2 - 1 



1/2 



2 v'l + 47™^ + 1 



= (^p^ - ^ ^; j . (71) 

When — 1 < Arn^j < 0, both frequencies are real. When Ani^j < — 1 both fre- 



quencies have an imaginary part. In the case 7 > 0, both terms ^yl + 4777?,^ ± IJ , 
are positive so that ujp is always real, while flp turns out to be imaginary when 

< + A'yrn? + 1^ /2'y. Those sectors of the theory with complex frequencies 
give rise to runaway solutions. 

From the canonical brackets 

{ip{x,t) ,'K^{x',t)} = iS{x - x') , {ip{x,t) jTT^ {x',t)} = i5 {x - x') , (72) 
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together with the foUowing equations of motion arising from (I66p 

(p^tp, = -TT^, vr^p = (1 + 27A) V' - TT^, , (73) 
7 

which allow us to express the remaining fields ip, tt^ and n^p in terms of Up, a^, , bp, 6^, , 
we get 



1 



-5' (P - P') , 



{ap,ap'} = |a]„a]^,| = {bp,bp,} = |?^]„6p.| = 0, 



(P-P') 



(74) 
(75) 
(76) 



which is the same result obtained in Ref. 1I6I with 7 — >■ —7 and a different definition 



for the auxiliary field. Let us recall that any Poisson bracket between ap,al, and 



bp, bp, is zero. With the above normalization the number operator is 



N = v^l + 47^2 J (fp [alap + b^] . 
In terms of the momentum space fields the Hamiltonian becomes 
H=^ d^p {{apal + alap) ujp + {bpbl + blbp) Vlp) . 



(77) 



(78) 



From here the canonical quantization is straightforward. The algebra of the creation 
and annihilation operators is 

1 



\ap, ap' 



1 



'5-'(p-p') 



-J^ (p-p') 



^1 + 47m2 



«P'"p' 



= 



and the normal ordered Hamiltonian operator results 

Hn = j d^P {alapOjp + blbpflp) , 
which acts on one-particle states according to 

iJ^4|0,0) = ^=^^4|0,0) 



HNbl\0,0) 



^1 + 47m' 



4|0,0) 



(79) 
(80) 
(81) 

(82) 

(83) 
(84) 



Due to the minus sign in Eqs. (|80l84p . the requirement of energy positivity implies 
a negative metric for states containing type b particles, spoiling physical unitarity 
and a consistent probabilistic interpretation. A normalized state with n particles 
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(s particles of type a, with momenta pi, ■■■■,Ps and {n — s) particles of type b with 
momenta Ps+i, ■■■,Pn) in the a = case is given by 

\'fa,biPi, P2,:;Pn)) = {i + 4^mY (n^iO (nr=.+i<.) |0,0). (85) 

The expression (j84l) shows that the minus sign in the commutation relations 
((80)) leads to negative contributions for the energy from the excitations of type b. 



4.2. The perturbative Hamiltonian constraint approach 

Here we go back to the self-interacting case a ^ 0. The expression for tt^, equation 
(|58p . is non homogenous in 7, and hence generates a Hamiltonian perturbative 
constraint of order 7" 

7"^^ = 0, (86) 

where n is the order of the perturbative approximation we want to achieve. Thus, 
the primary Hamiltonian density for a perturbative approach of order n in 7 is 

n = - (1 + 27A) V + (m^ - A - 7A2 + ^ 

+ tp-K^ — tta (m^ — A — 7A^ + aip^) ip + vj"'n^ (87) 

and the full set of constraints, primary and secondary, are given by (|59p and 
together with the chain 

$0 = 7"^^ - (88) 

*o - 7" (V- - vr^) - (89) 

$1 =7"-i(7(TO2_A + a(^2)(^ + 7rv;) ~0 (90) 

^5>^ 7«-i [(7 {A + m^ + 3aip^) + l) t(; - tt^] ~ (91) 



which are obtained by requiring 

{i/^ vj/"-™} = (J)"-™-!. (92) 

As a simple illustration, let us consider the case n = 1. The above chain produces 
the additional perturbative constraints 

$1 = 7 (m^ - A + a(y3^) (y3 + TT^ = (93) 
*i = (7 (A + TO^ + 3a^2^ + 1) ^ _ ^ (94) 

which have to be considered together with the original ones ([55]) and (|60p . 

The set ([55]) . (|60p . and ([M)) corresponds to four second class constraints, 
leading to the Dirac brackets 

{V,'^v}d = W'^v} - {</':*l}C'^^j^<I.i {^l.T^vl - {¥',X2}C^2*1 {'^'l''^¥'} 
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where Cab is the standard matrix of the Poisson brackets among the second class 
constraints. The final result in the [lp^t:^) reduced phase space is 

{^,^^}^ = l-7(m2-A + 3V-)^ (gg) 

together with the Hamiltonian density 

n ^vr^ (1 - 27A) TT^ + (1 + 7m2) <^2 _ 1 _^ 277772) ipI^Lp (l + 47777^) 

+ I + ^7aV® - la {iif^ (V(p)^ + 2(^3A(^) . (97) 

The final Dirac bracket (|96l) is not canonical. To express the Hamiltonian density 
in terms of canonical variables, maintaining the original field Lp, it is necessary to 
apply a non canonical transformation 

= (p, (98) 
(1+7(7772- A + 3a(^2^)7r^, (99) 



such that 



{^,^4^ = 1 + 0(7'). (100) 

In terms of this new momentum, which we call again tTj^ in an abuse of notation, 
the Hamiltonian density ((97|) becomes 



n = ivr,^ (1 - 27 (7772 + iaip^)) TT^ + ^777^ (1 + 7777^) (^^ 

- i (1 + 27 (7772 + atp^)) tpAip + 7 (1 + 47777^) + \-ic?^P. (101) 

It is straightforward to verify that this Hamiltonian density yields the same equation 
of motion ([56)) as the original Lagrangian density (|54)) . to first order in 7. To close 
the discussion, we obtain from Eq. (jlOll) the effective Lagrangian density to first 
order in 7 in the configuration space ip 

= T^ip^ - "H = - + I'm^ {p^f + 771^ (1 - 7771^) (p^ 

+ ^ (1 + 27777^) ip"^ + 2-iatp^ULp + a^-fip^^ . (102) 



This Lagrangian density has to be compared with the exact one, given by Eq. ([54| . 
The effect of the HOTD term in this effective Lagrangian has been to produce mod- 
ifications in the mass of the field ip, in the coupling constant of the self-interaction 
and has also generated two new interaction terms, one of order ^ip^ with a 
derivative coupling and another of order 7(^5^. Normalizing the kinetic term via the 
substitution 

(/? ^ = (1 + 7777^) ^, (103) 

we get 

Cd = - (1 - im") ^2 - _ ^„^3g^ _ ia'7^'- (104) 
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Before closing this subsection it is interesting to compare the final Lagrangian 
density (jl04p . obtained via the perturbative Hamiltonian constraint method, with 
different alternatives previously proposed to reduce (I54[) to an effective first order 
form. 

One such alternative is the double zero methotJ^. In this case we can use de 



equation for motion (j55p to generate such a term, which leads to 

Cdz = -^'pOip - ^(pm'^ifi + '^ifU'^f - jf'^ - ^7 {dip + m^ip + aip^Y , (105) 

and thus the effective Lagrangian to first order in 7 results identical to (jl02l) . so 
that after the redefinition (|103|) reproduces (|104p . 

Another possibility is to implement an appropriate derivative transformation in 
()54|) , to first order in 7, given by 

p^p = (^ip + ^iD^^ , (106) 

which yields the following first-order Lagrangian density to order 7 

£ - - i (1 + 7m2) ipD^ - im2^2 _ _ ^^a^3^^, (107) 

We can rewrite this expression in terms of the original field ip perturbatively, to 
first order in 7, by using equation (|55p 

0- (^V-^^l^V^ ^V+^lim'^f + aip^) ■ (108) 
In this way we get the following Lagrangian density, to first order in 7 
£g = — ^- + 7m^^ {fOip + (1 — 7m^) ip"^ 

+ ^ (1 + 2jm^) ip^ + 2jaip^nip + a'^jip'^'^ . (109) 

This Lagrangian density is identical to (|102p which also reproduces the one given 
by the double zero method. Summarizing, to first order in 7 we have 

Cd = Cdz = Cg- (HO) 



5. The two-particle scattering in the Bernard-Duncan scalar field 
with a interaction 

In this section we discuss the dispersion of two scalar particles, with a dynamics 
described by the Lagrangian density (|54p . which contains the interaction term 

C^nt = -^V^. (Ill) 

To test the perturbative Hamiltonian constraint approach we analyze this quantum 
process following two different approaches. In the first place we directly compute 
the scattering amplitude at first order in 7 using the results of the perturbative 
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Hamiltonian constraint method. After this we obtain the expression for the scat- 
tering amplitude exact in 7 from the Lagrangian (|54p . from which we derive the 
corresponding amphtude at first order in 7. We assume that the quantum Hamiho- 
nians are defined with the normal order product, so that the tadpole diagrams are 
not considered. 

It is interesting to study the dispersion of two scalar particles at first order in 7, 
which in the effective theory involves not only corrections to the original vertices, 
but also the new derivative vertex. To proceed with the first calculation we consider 
this process using the Feynman rules derived from the effective Lagrangian (|102p . 
It corresponds to the effective Hamiltonian (jl01|) . obtained using the perturbative 
Hamiltonian constraints approach, for which we have the following algebra for the 
creation and annihilation operators 



such that the normalized two a-particle in and out states are 

1 



(112) 
(113) 

(1- 277712) ■ ^^^^^ 
According to this, the first order contributions to the scattering amplitude (see 



($o«t| = 



(l-27m2) PI 
, . 1 




i (1 - 27m2) 
p2 — rm? (1 — 777*2 




&ia^ (pf + pf + pI 



-pI) 




Fig. 1. Feynman rules for the effective scalar theory. Note that now we have two additional 
vertices, a ip^ derivative vertex and a ip^ vertex. 
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Fig. [1]), are the ones given by the vertex 



3i {2?:^ a 



2 ■y/^^W^W^M'^ 

and the derivative vertex 



(1) _ 6j_(27r)a7m2 



'ft 



such that the total first order contribution results 



(P'l +P2 - -Pi --P2) , 



(1) 3»(27rf a(l-47m^) ^,^ , , 
Sf/ = 5 [Pi +P2-P1- P2) ■ 



2^ujp^ujp^ujp'^ujp'^ 



(115) 



(116) 



(117) 



At second order we also have two contributions (see Fig. [2]). One of them corre- 
sponds to one loop with two a* type vertices 

Sf^ - _ 9 "M^-47m^) +p^)+F ip[ ^p,) + F ip', - p,)) 

xS\p[+p'2-pi~P2), (118) 
^ (119) 



F (p) = J d*k 



^2 _ ^2 (1 _ ^^2) _ ^^2 _ ^2 (1 _ ^^2) ' 




Perm. 




Perm. 




Perm. 




Perm. 



Fig. 2. Feynman diagrams contributing to the two-particle scattering in the effective scalar theory, 
at first order in 7 and one loop approximation. 
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and the other to one type vertex and one derivative vertex 



xS*{p[+p'2^Pi^P2), (120) 



(2m^ + k^ + {p- kf^ 



1 



G{p)^ d^k^- — 2 . (121) 

' /c2-m2(l-7m2) (p-fc)^-m2(l- 777i2) 



The complete expression for the scattering amphtude in this approximation is 



sy^ = ' ' ' -r {p[ + P2 -pi- P2) 

Z^UJp^UJp^UJp'^UJp;^ 

X ((^^ + 7G)(,,+,,) + (F + 7G)(^,_^^) + (F + 7G)(^,_^^)), (122) 

r I 1 - (2m^ + k^ + {p ^ kf) 

iF + .G),,, = y _ _ ^,„.) ^ ^ (123) 



Note that Ff, contains a quadratic divergence and a logarithmic divergence. One 
should regularize this quantity in order to renormalize the theory. However, the point 
we are interested in is comparing the results for the scattering amplitude using two 
different methods, computing in the effective Hamiltonian theory constructed using 
perturbative Hamiltonian constraints in the usual fashion and computing in the 
exact HOTD theory, but approximating the propagator and restricting the initial 
and final states. For this reason we do not make explicit here the renormalization 
issues. 

We now consider the same process, but in the framework of the exact Lagrangian 
(|54l) . To compare with the results from the perturbative Hamiltonian constraints 
approach, in terms of the Fourier decomposition (|69p . we consider only the asymp- 
totic states corresponding to physical, a-type, particles, not including negative norm 
ghost contributions in the initial and final configurations. Thus, according to (|85p . 
the in and out states are respectively 



= ^l + 4jm^alal |0) , (124) 
($ont| = (0|ap;ap^yrT4^. (125) 



Hence, after adding all contributions up to first order in 7 see Figs. ([3]|4|), the 
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scattering amplitude at first and second order in a are respectively 



Si (27r)''a(l -47777?) 



S ip'l+P2-Pl -P2) 



(126) 



4? = ^ (^) I I d^x, (0| ap,ap.^ T (r^^ {x,) : -.^^ (x^) : ) a^a^ |0) 



9a^ (1 - 47777^ 



where 



2^w^ 



(Pi +P2 -Pi --P2) 

X f F (pi +P2)+F {p[ -pi)+F - pi) 



F (p) = / d^'fc 



7fe4 + fc2 - m2 ^ (fc _ p)4 _^ (fc _ _ ^2 ■ 



(127) 



(128) 



7P'* + p'^ — vr? 




-Gia 



Fig. 3. Fcynman rules for the HOTD Bernard-Duncan theory. 




+ Perm. 




+ Perm. 



Fig. 4. Feynman diagrams contributing to the two-particle scattering in the HOTD Bernard- 
Duncan theory, to first order in 7 and one loop approximation. 
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At first order in 7 they reduce to pT7)l and (fT^ . witlr F {p) = (F + 7G)(p) +0(7^). 
Both computations, the one based on perturbative Hamiltonians constraints and 
the one obtained from a perturbative expansion of the propagator of the exact 
Lagrangian, where we have used a procedure similar to the one applied by Weinberg 
in Ref. [HI yield the same result at the order 7 considered. This is evidence of the 
soundness of the effective Hamiltonian construction proposed here. 

6. Final remarks 

In this article we have presented an alternative method for constructing a consistent 
effective Hamiltonian formalism for higher order Lagrangians, which gives a correct 
low energy approximation. It includes higher energy scale effects under the form of 
perturbative corrections in the framework of a second order theory, so that it is free 
from the pathological behavior characteristic of HOTD theories. This method re- 
sults from an application of the standard Dirac procedure to deal with constrained 
systems which, together with the constraints generated by the definition of the 
momenta and their consistency conditions, incorporates a new set of perturbative 
constraints generated by the inhomogeneous relations in the perturbative param- 
eter contained in the set of the original momenta definitions and constraints. The 
addition of these new constraints projects the dynamics of the original HOTD sys- 
tem into a stable subspace consistent with the chosen order of approximation. Our 
method is purely algebraic and does no require to solve any system of differential 
equations. Also, it can be directly applied to gauge systems without the necessity 
of initially fixing the gauge. Moreover, we conjecture that the brackets obtained 
from this procedure in phase space, when rewritten in the CV space via the corre- 
sponding equations of motion, would provide a solution for the system of differential 
equations (|A.10p . which is the starting point of the Eliezer-Woodard formulation 
131. The resulting second order effective Lagrangians produce unitary theories after 
quantization. The study of the possible breaking of Lorentz covariance induced by 
the perturbative description is beyond the scope of the present work 

At the classical level this Hamiltonian formalism yields canonical equations of 
motion equivalent to the usual perturbative approximation for the Lagrangian equa- 
tions of motion of the exact HOTD theory. Furthermore, from this canonical con- 
struction we can also derive a well behaved effective Lagrangian formalism. At the 
quantum level we recover the results obtained using the perturbative expression for 
the propagators, provided that the asymptotic state space is restricted to the phys- 
ical one, ruling out the ghosts. The difference with this last approach is that now 
we have a well defined canonical formalism, instead of a non-consistent one where 
ghosts must be forbidden because asymptotic states and approximate propagators 
are used. This new formulation provides a consistent effective theory, which allows 
the implementation of all the usual manipulations for the construction of a quantum 
field theory. 

An important characteristic of our approach is that it makes no use at all of the 
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Lagrangian equations of motion or any Lagrangian constraints to obtain the cor- 
responding effective theory. This is a significant difference with previous proposals, 
particularly the methods presented in Refs. l4land[T2l Moreover, our formalism can 
deal with singular theories in a straightforward way, since we only have to add the 
corresponding non perturbative relations to the total set of constraints and work 
within the Dirac framework in the usual way. Another characteristic is that the pro- 
cedure is not an iterative one. At the start we choose the order of the approximation, 
which is defined by the chosen perturbative constraints, and the Dirac algorithm 
converges directly to the corresponding effective theory. 

As working examples, we have applied the method to two paradigmatic models 
in HOTD theories: the Pais-Uhlenbeck oscillator and the Bernard-Duncan scalar 
field, in which we include a (p'^ interaction. The first example clearly shows the 
main features of the perturbative Hamiltonian constraint approach, and allows us 
to express the perturbative constraints in a very simple and closed way, in terms of 
Fibonacci polynomials. In the other example, closed expressions for the constraints 
are much more difficult to write and for this reason we restricted the construc- 
tion only to first order in the perturbative parameter. In both cases the canonical 
equations of motion are equivalent to the Lagrangian perturbative ones. 

From the Hamiltonian thus constructed we can obtain a second order effective 
Lagrangian. The quantum theory can be obtained either from the Hamiltonian for- 
malism or from the effective Lagrangian approach, via the Matthews' theorem which 
is well established for Lagrangians of the form here obtainecfl^. The calculation of 
the two-particle scattering shows that the effective theory constructed on the basis 
of perturbative Hamiltonian constraints gives the same results as the exact higher 
order theory, provided that in this last theory the space of states is restricted to 
the physical one, ruling out ghost states, and that the propagator is considered in 
terms of a perturbative expansion. Similar results are obtained in Ref. 1181 although 
there a direct substitution of leading order equations of motion was performed in 
the Lagrangian to get the effective theory. As it has been clearly shown in Ref. [231 
this procedure can be considered to be correct only when it effectively results in 
the addition of a double zero term to the original Lagrangian, or when a suitable 
equivalent derivative field transformation is found. 
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Appendix A. 

In this Appendix we explore the connection between Lagrangian and Hamiltonian 
perturbative constraints, in the context of a simple example of HOTD theory. Let 
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US consider the Lagrangian 

= ^ i<i^f - iltf + l^^iAi<i3- j = 1, 2 (A.l) 

For simplicity, we will construct the Hamiltonian formalism only to first order in 7. 
The exact equations of motion are 

q^+le^Jqf +Lo'^q^^Q (A.2) 

which, to first order in 7, reduce to 

iii = -fu'^eijCij - uj'^Qi (A. 3) 

In this case, the approach of Jaen, Llosa and Molina ^l^l and the one of Eliezer and 
WoodardSlare equivalent, as stated in Ref. El For this reason we will consider only 
the second one, which is more adequate to establish the comparison. Following Ref. 
IHwe calculate the Noether energy 

E{q, q) = ^qf + ^oj^qf + l^rM,, (A.4) 

which is taken as the time evolution generator in the coordinate- velocity (CV) space. 
In principle E{q, q) should be a projection on the CV space of the Hamiltonian which 
lives in a not yet known phase space. To first order in 7 the above equation yields 

El ]^ql + ^t^^gf + 7^^eij9j9j- (A.5) 

The method further assumes the existence of a fundamental bracket {qi , qj } in 
terms of which we can describe the temporal evolution, such that 

, , r 1 (^El , dE\ r K n\ 

q, = {qr,Ei} = {qi,qj} h {qi,qj} (A.6) 

dqj dq-j 

r • 1 r ■ 1 dEi dEi / A -T\ 

q^ = {qr,Ei} = {qt,qj} h {qi,qj} (A.7) 

dqj dqj 

These brackets must also satisfy the consistency conditions 

{q^: qj} + qj} = ^ ^j} , (A.8) 

and thus we obtain the following system of differential equations defining the basic 
objects in our example 



= {u'^qj + 7w^ejfegfc) {g,, q^} + (q^ - -^uj'^e.juqk) q^} , 
q^ = {uj'^qj + juj'^ejkqk) (^-^ {qi, qj} ~ 

+ {qj - J^^jkqk) + 2^') {9m 9,} , (A.IO) 
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where cji is given by Eq. (jA.3p . Following this approach, it is necessary to find the ad- 
equate solution to the above system, which in this case must be analytical in qi and 
Qi. Next, from this solution the relation between canonical momenta, coordinates 
and velocities must be inferred. Once this is achieved, the canonical formalism to 
first order in 7 is constructed by projecting the energy Ei together with the brack- 
ets among coordinates and velocities in the phase space just defined. Additional 
simplification of the system (jA.10[) can be achieved by making the following ansatz 

{qt,qj} = lAij, {qi-Aj} = St, + -fBij, (A.ii) 

to first order in 7. Even with the above simplification, some guess work has to be 
done in order to solve the Eqs. (|A.10p . It is clear that the complexity of the basic 
equations (|A.10|) will rapidly increase either when higher order approximations are 
considered or when more complicated systems are studied. We consider this as a 
shortcoming of the method proposed in Ref. IH 

For the purpose of comparing the method based on Lagrangian constraints, in 
the Eliezer and Woodard incarnation, with the one based on the Dirac approach, 
the expressions already obtained are enough. Now we will deal with the problem 
using the Dirac method augmented with perturbative Hamiltonian constraints. To 
apply the procedure we first rewrite the Lagrangian (jA.ip in first order form by 
introducing the coordinates Zi via the auxiliary coordinates Xi 

1 2 1 ? 2 1 

^9=2 ^ 2^ + 2'^i3^i^3 + i^i ~ 9*) ■ (A. 12) 

The canonical momenta are 

7 

Pq, = -Ai, PAi = 0, TT^^ = -£jiZj. (A. 13) 

So we have six exact primary constraints, plus two pertubative ones 

JTTz, = 0. (A.14) 

Following with the Dirac method we demand the consistency of the constraints 
under time evolution and we finally arrive to a Dirac Hamiltonian. Once the second 
class constraints are imposed as strong relations it reduces to 

He = \v\ + (A.15) 

with the Dirac brackets 

{q%-,Vj)B^ ^n-- te^PjIn^O, {qi-,qi)o^i^-io (A.16) 

To project this formalism in the CV space and compare with the Eliezer- Woodard 
approach, we perform the transformation ((/,(/) (q-p) which is given by the 
canonical equations of motion 



q^ = {qi,Hc} =Pt+ -fJ^ei-jq-j 



(A.17) 
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Projecting He we obtain 

He = \qi + \^'^qi + 7w^e»jft9j-, (A.18) 

which coincides with the Noether energy restricted to the C V subspace (jA.Sp . From 
the Dirac brackets (|A.16I) and the transformation (jA.17[) we can compute 

9j} = + l^'^'ijkqk) = k]-, (A.19) 

^j} = {Vi + Jc^'^e^mqni,P] + juj^ejkqk} = 2juj'^ei-j. (A. 20) 

It is straightforward to verify that this set of brackets is indeed a solution of the set 
of diferential equations (jA.lOp . 
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